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Abstract 

We study the time evolution of an N-component model of bicontinuous mi- 

croemulsions based on a time dependent Ginzburg-Landau equation quenched 

from an high temperature uncorrelated state to the low temperature phases. 

The behavior of the dynamical structure factor C(k,t) is obtained, in each 

phase, in the framework of the large-iV limit with both conserved (COP) 

and non conserved (NCOP) order parameter dynamics. At zero temperature 

the system shows multiscaling in the unstructured region up to the tricriti- 

cal point for the COP whereas ordinary scaling is obeyed for NCOP. In the 

structured phase, instead, the conservation law is found to be irrelevant and 

the form C(k,t) ~ t a / z f{(\k — A^i 1 / 2 ), with a = 1 and z = 2, is obtained in 

every case. Simple scaling relations are also derived for the structure factor 

as a function of the final temperature of the thermal bath. 
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Self-assembled fluids are mixtures of water, oil and surfactants displaying a richer ther- 
modynamic behaviour than ordinary fluids because they appear to have a tendency to order 
and form patterns. In fact, if the surface tension can be made sufficiently low, by increasing 
the surfactant concentration, one observes a microemulsion phase stabilized by spontaneous 
curvature and entropy of mixing, whereas at higher values of the surface tension one obtains 
oil-rich or water-rich regions. By increasing further the surfactant concentration, however, 
the fluid acquires an internal order as revealed by the several structured phases. 

So far most of the experimental and theoretical studies have been devoted to the in- 
vestigation of the equilibrium properties of such systems whereas the properties far from 
equilibrium of complex fluids are not so well understood [0J. In the present Letter, by gen- 
eralizing to a large number of components the GL model for oil-water-surfactant mixtures 
put forward by Gompper- Schick ||, ||, we study the dynamics after a quench. 

By associating an order parameter with the difference in the density of water and oil 
molecules we build the GL functional as the sum of a local term plus a non local interaction, 
which is mathematically convenient to represent as a gradient expansion. The model is 
described by the following GL hamiltonian: 

H[<t>{x)\ = J dx[\{y^f + b -(v<t>) 2 + ^(</>) 2 W) 2 + r -(<t>) 2 + ^((</>) 2 ) 2 ] (i) 

where r and g (with g > 0) are the quadratic and quartic couplings of the GL theory. The 
term containing the squared Laplacian represents a curvature energy contribution, whereas 
the second and the third term are proportional to the surface energy. The peculiarity of the 
self-assembled fluids is reflected by the presence of higher order derivatives in the gradient 
expansion and by the sign of the coefficient b of the square gradient term, which can be 
either positive as in ordinary fluids or negative according to the value of the thermodynamic 
control parameters, as found from the fits of the experimental data |3J, while the constant 
w is assumed to be positive. 

Incidentally, our results turn out to be of interest for a classical problem of relaxation 
dynamics, apparently quite different as the hydrodynamic fluctuation at the convective 
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instability, the so called Swift-Hohenberg model 0- 0. 

Within the large N limit the equation of motion for the order parameter and its fluctu- 
ations can be written in a closed form and we can extract explicitly the exact asymptotic 



behaviour ||- JTT]|. A continuum description of the model begins with associating the mod- 
ulus of the vector order parameter <p(x, t) = (<fri(x, t), .., 0at(x, £)) with the difference in the 
density of water and oil molecules and assuming that its evolution towards equilibrium is 
described by the Langevin equation: 

°' ( X - ' ' -T(x)—^-—H[cj>(x)} + Va (x, t) (2) 



dt 5(f) a (x,t) 
Here ff(x,t) represents a Gaussian white noise with zero average and < i] a (x,t)f]f3(x' ,t) >= 
2TfT(x)5 a ^5(x — x')5(t — t'), where Tf is the temperature of the final equilibrium state and 
the kinetic coefficient takes a constant value T for non conserved order parameter (NCOP), 
whereas it is given by —TV 2 in the conserved case (COP). 

In order to study the properties of the system one considers the long-range behaviour of 
the equal time real space connected correlation function C(r,t) =< (fi a (R + r,t)(f> a (R,t) > 
and its Fourier transform, the structure factor, C(k,t), which are both independent of the 
index a due to the internal symmetry. 

In the large iV-limit from eq. (fj) one derives the following evolution for the structure 
function : 

j-C{k, t) = -2Tk p [k 4 + B{t)k 2 + Q{t))C{k, t) + 2k p TTf (3) 

with p = for NCOP dynamics and p = 2 for COP. We have introduced the two auxiliary 
functions B(t) and Q(t), which have to be determined self-consistently through: 

B(t) = b + wS (t) (4) 

and 

Q(t) = r + gS (t) + wS 2 (t) (5) 

where the integrals 



s °^ = L (6) 

J|fc|<A {2TT) a 



and 

d d k 



= I K v^r d k 2 C(k,t) (7) 

J\k\<A [2ll) a 

contain a phenomenological momentum cutoff A. 

At equilibrium ( see also ||12||), when r < 0, the system displays a low temperature or- 



dered, ' magnetic' phase with non vanishing order parameter < cb > and a high temperature 
disordered (paramagnetic) phase as shown in Fig. [l]. In order to draw the phase diagram of 
the system we consider the equilibrium value of the structure factor C eq {k) which reads: 

T 

where we have introduced b r = limt_ >00 (6 + wSo(t)) and D = limi__ )00 (r + gSo(t) + tuS^t)). 
The 'ordered phase' < cb ^ > is bounded from above by the line of 'magnetic' critical 
points T c (b) shown in Fig. [I] . It terminates at the tricritical point (pL = rw/g, T tr ), 
located at T tr = in d < 4. Points below such a curve correspond to a vanishing value of 
the parameter D. In other words the line T c (b) separates a high temperature phase with 
finite correlation length and finite fluctuations from a low temperature magnetic phase with 
infinite correlation length and divergent fluctuations at zero wavevector, due to the presence 
of massless Goldstone modes. 

Within the disordered phase the two-particle correlation function displays two different 
behaviors according to the sign of the discriminant A = [D — h-]. In fact, the vanishing 
of A identifies the so called disorder line, i.e. the borderline between a regime (since it is 
not a proper thermodynamic phase) with monotonically decaying correlations, which is the 
analogue of a paramagnetic phase, and a phase with oscillatory-decaying correlations, which 
in the present model corresponds to a microemulsion phase. The locus where the structure 
factor, C eq (k) starts developing a peak at finite wavevector k, when b r (Tf) = 0, is named 
Lifshitz line. It is contained in the microemulsion phase and terminates at the tricritical 



point. To summarize, for positive values of A, C(r) ~ 6 r - - sin(k m r) with | = [^\^D + y] 1 / 2 
and with k m = [^y/~D — ^] 1//2 . For negative values of A, the correlation decays monotonically, 
showing an Ornstein-Zernike behavior. Finally along the semiaxes Tf = and b < rw/g, 
A = and the structure factor diverges at a finite wavevector. Such a line is the remnant 
of the lamellar phase, which is unstable for any finite temperature, because the system does 
not support topological defects in any dimension d less than N. 

We shall consider now the phase ordering of the present system following a quench from 
an uncorrelated high temperature state for different choices of b r . In order to characterize the 
quench process one needs to specify the the set of parameters b, w, r, g and initial conditions 
for the field; this can be done by assigning the structure factor at t — 0. One usually assumes 
the initial state to be uncorrelated and chooses the initial condition C(k, 0) = A = const. 
In the following we shall consider r < 0. The asymptotic form of C(k,t), for the Tf = 
case, can be easily calculated by imposing the matching of S (t) and 5 2 (i), calculated by 
means of the general solution of eq. (|j) with their equilibrium constant values. 

In the NCOP case and for positive values of b r , since the curvature term is asymptotically 
irrelevant, the structure factor has the usual scaling form as for simple fluids C(k,t) ~ 
t a l z f(kt l l z ), with f(x) = e~ x2 , a — d and z — 2, for a quench below the critical temperature 
(Tf < T c ). The exponent z = 2 is, in fact, related to the growth of the domain size 
L(t) = (2rt) 1 / 2 , which is controlled by surface tension. In the COP case one observes, 
instead, the multiscaling behaviour of the structure factor |J, due to the existence of two 
marginally different scaling lengths, L(t) = (2rt) 1 / 4 and k^(t) ~ (t/ logt) 1 / 4 , where k m (t) 
is the position of the peak of the structure factor C(k, t). This results in a multiscaling form 
for C(k,t) ~ (k 2 ~ d L 2 )^ k l km) with <j>{x) = 1 - (1 - x 2 ) 2 . 

Of particular interest is the study of a quench at the tricritical point, where the 'mag- 
netic', the 'paramagnetic' and the 'lamellar' phase meet. This occurs when b r vanishes, at 
Tf = 0, i.e. in correspondence of the bare value b^ = rw/g. 

We study, first, the growth of the structure factor upon approaching Tf = along a 
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particular Lifshitz line, by considering the simpler case w = 0, which yields bi = 0. For 
the NCOP at T f = the structure factor reads C(k,t) ~ A e~ 2rfc V/ 4 , i.e. a scaling form 
characterized by a = d, z = 4 and scaling function fix) = e~ x ' 4 . The typical domain size 
evolves according to Lit) = (2rt) 1//4 , being controlled by the curvature rigidity and not by 
surface tension. For the same choice of parameters in the COP case, instead , the growth 
process depends on two distinct lengths Lit) = (2rt) 1//6 and k^(t), where k m (t) is the 
wave vector where C(k,t) reaches its maximum and varies in time as k m (t) ~ ^d]M^i/6_ 
The same mechanism which determines the multiscaling behaviour of the structure factor 
upon approaching an ordinary critical point leads to the multiscaling form for C(k,t) ~ 
(A4- d L 3 )^ fe / fe -) with 4>{x) = (3x 2 - x 6 )/2. 

For the choice b = 0, w = and Tf > we have investigated numerically the behaviour 
of the solutions of eq.([jj), with NCOP dynamics. As shown in Fig. || by plotting C(k m ,t) 
times Tj^ 4 d ^ against time multiplied by d ^ one obtains a remarkable data collapse 

over several decades. This suggests the existence of a typical relaxation time r which scales 
with the temperature Tf according to the law r ~ TJ 4 /( 4-a! ) . Such a dynamic scaling relation 
can be deduced heuristically, by matching the typical value of the maximum of the structure 
factor, which grows as ~ t d l 4 up to the characteristic time r, with the equilibrium value 
Ceqikm) at Tf which is proportional to Tf d ^ 4 d \ 

Let us consider a slightly different case, where b = rw/g, i.e. a quench to the tricritical 
point at (6 ^ 0, Tf = 0). Here the two lengths A;" 1 and £, discussed above, diverge upon 
approaching equilibrium. One observes a crossover from an initial regime characterized by 
a structure factor with a peak at finite values of k to the true asymptotic regime displaying 
the ordinary Bragg peak at k = at late times. For the NCOP dynamics one finds from 
the self-consistency conditions eqs. (Q) and ([5]) that asymptotically 2TB (t) ~ at 1 / 2 , with 
a > 0, and 2TQ{t) = — G?log(i)/4. Therefore the structure factor assumes at late times 
the scaling form: C(k,t) ~ t d / 4 exp(— x 4 + 2ax 2 ) where x = (2rt) 1//4 fc. Notice the unusual 
feature, for NCOP, of the maximum of the structure factor at finite wavevector moving 
slowly at the rate ~ t -1 / 4 towards the origin . In other words initially B(t) < (for times 
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such that So(t) < \b\/w), and the system promotes fluctuations with k m (t) > 0, which are 
eventually suppressed by the growing of correlations which cause k m — > 0. However, the 
two characteristic lengths L(t) and k^{t) vanish at the same rate and standard scaling is 
observed. 

A novel behavior is observed when one crosses the Lifshitz line and b r < since asymp- 
totically the structure factor C(k,t) displays a peak at finite wavevector, which is given by 
k^it) = —B(t)/2t . The main difference with the cases treated above is the finite value of 
fc" 1 in the late regime. 

For the NCOP dynamics at Tf = a saddle point estimate of the integrals yields: 



S (t) ~ K.e-^e^l^^f-^i-Bit)]- 1 / 2 
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(9) 



S 2 (t) ~ K d e- 2T ^e^[-^]^l 2 [-B{t)]- l l 2 (10) 

with K d = [2 d - 1 n d / 2 T(d/2)}- 1 . Upon requiring that So(t) and S 2 (t) reach a constant value 
at equilibrium we obtain: 

Q(t) ^__L log[ _ B(t)] + (^ log[ zM, 

For consistency, since in this region for large times B{t) < and — > k 2 m we find 
asymptotically B(t) ~ —2k 2 m t. After the quench from a high temperature initial state 
the structure factor develops a peak at finite values of k, corresponding to a quasi-ordered 
pattern, with periodicity given by -jp. The height of the peak grows in time according to 
the power law t 1 ! 2 . This is due to the competition between curvature and surface tension 
controlled fluctuations which eventually lead to a layered structure at Tf = 0. Within this 
regime we find the following asymptotic form of the structure factor 

C(M) ~ v / texp-2r(A; 2 + ^) 2 t (12) 

which obeys the scaling form C(k,t) ~ t a l z f(\k — k m \t x l z ') with a — 1 and z — 2. This 
unusual lack of dependence of the structure factor on the dimensionality, d, reflects the fact 



that the largest fluctuations characterizing the ordering process have a finite wavenumber 
and thus the density of state contribution does not manifest itself. We have considered also 
the finite temperature behavior of eq.(|3|) and found data collapse for the curves representing 
the height of the peak times Tj versus the time multiplied by T 2 . In fact, the height of the 
peak grows according to eq.(|T2"D as t 1//2 , while for Tf > a simple calculation shows that the 
structure factor settles at a finite value C eq (k m ) ~ Tj . By matching the two curves one 
finds a relaxation time r proportional to Tj 2 . 

In the case of the COP dynamics the only qualitative difference observed is the weak 
dependence of the peak position on time, since k m (t) approaches as log(t)/t its equilibrium 
value. In this case, therefore, the conservation law turns out to be an irrelevant constraint 
since it acts at k = 0, while the most important fluctuations are grown at finite k. 

A natural question is to ask ourselves how generic is the t 1 ^ 2 behaviour we have observed, 
namely whether the present results for iV — > oo have any relevance with respect to other 
models. We believe that our findings are consistent with the very late time behavior (z = 2) 
observed in numerical simulations by Elder et al. []7| in a study of the NCOP Swift-Hohenberg 
equation with N = 1 and d = 2 within the structured phase. On the other hand, a part 
for a special choice of parameters, we do not find any evidence in our model of their growth 
exponent z = 4 related to curvature relaxation, a feature probably due to the lack of 
topological defects. We finally remark that the dynamical equations we solved are very 
closely related to those derived by employing the dynamical Hartree approximation for 
N=l Q. 
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FIGURES 

FIG. 1. Phase diagram for r < and d = 3 

FIG. 2. Data collapse for the peak of the structure (T f = 0.1, 0.05, 0.01, 0.005, 0.001) factor 
at different temperatures for d = 3, b = 0, r = —0.1 , g = 1 , w = 
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